Nonlinear localized modes in two-dimensional electrical lattices 
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We report the observation of spontaneous localization of energy in two spatial dimensions in 
the context of nonlinear electrical lattices. Both stationary and traveling self-localized modes were 
generated experimentally and theoretically in a family of two-dimensional square, as well as hon- 
eycomb lattices composed of 6 x 6 elements. Specifically, we find regions in driver voltage and 
frequency where stationary discrete breathers, also known as intrinsic localized modes (ILM), exist 
and are stable due to the interplay of damping and spatially homogeneous driving. By introduc- 
ing additional capacitors into the unit cell, these lattices can controllably induce traveling discrete 
breathers. When more than one such ILMs are experimentally generated in the lattice, the interplay 
of nonlinearity, discreteness and wave interactions generate a complex dynamics wherein the ILMs 
attempt to maintain a minimum distance between one another. Numerical simulations show good 
agreement with experimental results, and confirm that these phenomena qualitatively carry over to 
larger lattice sizes. 

PACS numbers: 05.45.Yv, 63.20.Pw, 63.20.Ry 



Introduction. It has long been known that solitons 
emerge as classes of solutions to many nonlinear (lattice 
and partial) differential equations described chiefly by 
one spatial dimension; prominent examples are the KdV 
equation, the sine-Gordon or the nonlinear Schrodinger 
equations. In two dimensions, quasi one-dimensional lo- 
calization patterns can often still occur but robust 
two-dimensional (2D) localization in continuous media is 
rather atypical (see Ref. and references therein). How- 
ever, it is well-known that discreteness of the underlying 
medium can help to stabilize such localized solutions even 
in higher dimensions Alternatively, one can exter- 

nally enforce or introduce a periodicity in the form of a 
regular modulation in some property of the continuous 
(e.g. in optical photorefractive or atomic Bose-Einstein 
condensate) media @, Q , thus again breaking continuous 
translational symmetry. 

Here we show that two-dimensional discrete breathers, 
also known as intrinsic localized modes (ILMs), exper- 
imentally exist and are stable in the context of two- 
dimensional, damped-driven electrical lattices. We char- 
acterize these breather and multi-breather states in pa- 
rameter space, and we compare to numerical simulations 
and stability analysis. Finally, we focus particularly on 
versions of these lattices that support moving breathers. 
Discrete breathers have been considered in a variety of 
other settings experimentally including (but not lim- 
ited to) micromechanical cantilever arrays [7| , Josephson- 
junction ladders (JJLs) @], granular crystals of beads in- 
teracting through Hertzian contacts [9(, layered antifcr- 
romagnetic crystals [Io| . halide-bridged transition metal 
complexes ll| , and dynamical models of the DNA double 
strand fl2| . Yet, in most of these examples, the coherent 



structures are effectively one-dimensional. Even when 
higher dimensional (as is e.g. possible in optical waveg- 
uide arrays or photorefractive crystals Q), the states 
are typically stationary. Hence, the experimental abil- 
ity to systematically generate discrete breathers in two- 
dimensional electrical lattices and, perhaps especially, to 
control their mobility launches a new dimension in the 
modeling, and the theoretical and experimental under- 
standing of such states. This may be of broader interest 
in other areas as well, such as JJLs (see e.g. the theoret- 
ical proposal of Ref. |13j ) and two-dimensional granular 
crystals (see e.g. the recent experimental realization of 
Ref. [3l)- Another attractive feature of our lattices is 
their potential square or honeycomb geometry, especially 
since the latter has been a point of intense investigation 
at both the linear and the nonlinear level. This is due 
to its conical diffraction and Dirac (diabolical) points ex- 
amined intensely in both the physical (l5l . [Toj ] and math- 
ematical [TtJ communities, as well as due to potential 
connections with graphene nano-ribbons; see e.g. (l8l.ll9j 
and references therein. 

Experimental and theoretical setup. The experimental 
system under investigation is a class of two-dimensional 
electrical lattices of either honeycomb or square geome- 
try, as shown in Fig [TJ These electrical lattices can be 
considered as a set of single cells representing nonlinear 
LC oscillators, each composed of a varactor diode (NTE 
618) with a nonlinear capacitance C(V) and an induc- 
tor L 2 = 330 fiR. These single cells are then coupled 
at point V by inductors L\ = 680 /iH and are driven 
by a single sinusoidal voltage source V(t) via a resistor 
R = 10 kO, with amplitude Vd and frequency /. We 
study two different unit cell versions, with the only dif- 
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FIG. 1: Left: Schematic circuit diagrams of the basic geom- 
etry of the two electrical transmission lines (honeycomb and 
square), where black points represent single elements (right), 
with or without the block capacitor Cf. 



ference being the presence of a block capacitor between 
the diode and the coupling inductor in one of them, as 
explained in Ref. [2(J. The effect of the block capaci- 
tor is to make ILMs mobile; here we study two different 
values of block capacitances, Cf = 1 /iF and C/ = 15 
nF. The experimental lattices were comprised of 36 ele- 
ments with periodic boundary conditions, and node volt- 
ages were measured at a rate of 2.5 MHz using a 
multichannel analog-to-digital converter. In the linear 
limit of small-amplitude plane-waves, we have obtained 
the dispersion relation, which yields as the lowest fre- 
quency (uniform mode) f m - m = 312 kHz, and the high- 
est frequency / max — 689 kHz in the square geometry 
configuration, and / max — 617 kHz in the honeycomb 
lattice. Notice that this highest-frequency mode in the 
linear spectrum is now above the second-harmonic of the 
uniform mode in square geometry lattice, and below it, 
but very close, in the honeycomb configuration. 

Using basic circuit theory, in the block capacitor case, 
the dynamics of the lattice can be described by the equa- 
tions [U, 
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where the sum is taken over first, q, neighbors, q = 4 in 
square lattice case, and q = 3 in honeycomb one. Also, 
the following dimcnsionlcss variables were used: t = LUot: 
in,m = {Iv — I2) / {woCoVd) , where is the full current 
through the unit cell, and I2 the current through the in- 
ductor L2, both corresponding to cell (n, m); v = V/Vd, 
used with superscripts A and B; V 7 f m is the measured 
voltage at node (n,m); V r ^ m is the voltage at an inter- 
mediate point between the varactor and the block capac- 
itor; fl = 2irf/oj Q , ujq = l/y/L 2 C ; i D = I D / (ui C V d ) , 




FIG. 2: (color online) Comparison between the experimental 
((a) and (c)) and numerical ((b) and (d)) profiles of station- 
ary one peak breathers in a honeycomb lattice, (a) and (b), 
and a square lattice, (c) and (d). The insets show the Flo- 
quet numerical linearization spectrum, confirming (since all 
multipliers are inside the unit circle) the stability of solu- 
tions. The driver has amplitude and frequency Vd = 1.5 V 
and / = 283 kHz in the honeycomb lattice case, and Vd = 2 
V and / = 290 kHz (experimental) and / = 301 kHz (numer- 
ical) in the square lattice case. 



where Id is the current through the varactor diode; 
c = C(V)/C , where C = C(0), and C(V) is the nonlin- 
ear capacitance of the diode. Also, a phenomenological 
dissipation resistor, i?/, was included in the model to bet- 
ter approximate the experimental dynamics, and its value 
has been determined by matching the numerical and ex- 
perimental nonlinear resonance curves corresponding to 
a single element [2l| . A simplified set of equations corre- 
sponding to the non-block capacitor case can be obtained 
by considering the limit Cf — > 00. In all cases, the ratio 
L2/L1 characterizes the strength of the "effective" dis- 
creteness of the system. 

Results on standing/traveling discrete breathers in 
square/honeycomb 2D lattices. To generate discrete 
breathers, we have employed the well-known modula- 
tional instability of the driven uniform mode Q . In the 
case without block-capacitors, we obtain stationary and 
stable one-peak breathers, as shown in Fig. [2J These lo- 
calized modes arc robust, persisting as long as the driver 
remains on. 



Numerics are generally found to be in good agreement 
with experimental results. Yet, it can be seen in Fig. [2] 
that the amplitude of the one-peak breather is slightly 
higher in the numerics than in the experimental data. 
We attribute this slight mismatch to energy dissipation 
in the coupling inductors; these are stronger in the square 
lattice setting than in the honeycomb geometry, and they 
arc not taken into account in our simplified model of 
Eq. (jTJ) . When this dissipation is included in our model, 
by means of some phenomenological small resistances in 
series with inductors L\, this further improves the agree- 
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FIG. 3: (color online) Subharmonic breather corresponding 
to a square lattice, Vd = 6.4 V and / = 600 kHz. (a) Experi- 
mental spatial profile, and (b) experimental and (c) numerical 
temporal oscillations of three different nodes: the ILM center 
(black), the first neighbor (blue), and the node furthest away 
from the ILM center (red). Note that the oscillations of the 
node furthest away from the ILM center display local maxima 
corresponding to the natural period T of the ILM at twice the 
driving period T d — 1/f. 



ment for the ILM amplitude. For even higher amplitudes 
in the experiment, multi-peak breathers have been ob- 
served, which can also be captured by the model. 

As in the one-dimensional line [22j . we again ob- 
serve subharmonic breathers (experimentally and theo- 
retically) in two dimensions. In Fig. |31 the peak of the 
breather oscillates with a frequency /ilm = //2. No- 
tably, in two dimensions it is generally more difficult to 
stabilize breathers via subharmonic driving — the driv- 
ing conditions have to fall into a fairly limited region 
of parameter-space. Further illustrating the fragility 
of subharmonic breathers in two dimensions, we note 
that the introduction of block capacitor destroys sub- 
harmonic breathers — a qualitative difference with the 
one-dimensional lattice for which mobile subharmonic 
breathers were reported. The fact that the frequency of 
the subharmonic driver, namely twice the ILM response 
frequency, is now positioned within, or very close to, the 
linear dispersion band of the 2D lattice may explain this 
fragility. Our numerical simulations indicate (results not 
shown here) that in larger lattices, multi-peak subhar- 
monic breathers exist, and such subharmonic breathers 
appear to be more robust in the absence of block capac- 
itors. 

When a block capacitor of 1 fiF is placed in scries with 
the diode in the honeycomb and square lattice, the static 
breathers become mobile, as was the case for the one- 
dimensional chain [2(J. Nevertheless, and contrary to the 
simpler situation in one dimension, where a clear direc- 
tion of movement arises, the motion through the lattice 
appears to be more complex, as shown in Fig. 0] In gen- 
eral, numerical simulations obtained solving Eq. ([1} show 
fairly good agreement with experiments and yield infor- 
mation about the intrinsic nature of this complex mo- 
tion, where a strong sensitivity to small inhomogeneities 
is seen. The transition from one node to the next is 
such that where after a number of periods a noticeable 
asymmetry develops in the ILM profile; this leads to 
two neighboring nodes attaining equal amplitude, and 



finally the ILM becoming centered (initially between two 
sites and then) on the next node. This sequence is now 
also observed in two dimensions. Notice that the mo- 
bility problem has been argued to be quite important 
in other non-square lattices, such as the hexagonal one 
where breather mobility was proposed as responsible for 
the presence of dark lines in natural crystals of muscovite 
mica 1231. a nd reconstructive transformations in layered 
silicates [24| . In our square electrical lattice, we generally 
observe (cf. Fig. 2]) a directed motion which is interrupted 
by (longer) intervals of localization as the wave struggles 
to overcome the well-known Peierls-Nabarro (PN) bar- 
rier Q. Note that given the small size of the lattice, in 
addition to the role of inhomogeneities, small amplitude 
residual excitations ("phonons") are also important in 
directing the breather motion. It is also worth mention- 
ing that, although the movement along the square lattice 
usually happens along the edges of the lattice, some tran- 
sitions are also observed to happen, both in the experi- 
ment and the numerics, along the diagonals that have a 
slightly higher PN barrier. Despite the complex nature 
of the motion of the breathers, the organizing principle 
obeyed in both one and two dimensions is that the ILM 
never hops back to the node it occupied prior to its cur- 
rent location due to some residual capacitor-charge im- 
purity associated with that node. In one dimension, this 
principle necessarily gives rise to orderly, uni-directional 
motion. In two-dimensional hexagonal lattices there are 
still two choices available to the ILM, and in a square lat- 
tice there are three. Thus, the motion does not have to 
be uni-directional. It would be interesting to statistically 
measure the "diffusion" associated with the apparent ir- 
regular trajectories displayed by the moving breathers. 
Although the system is deterministic, due to small de- 
fects and residual excitations, the breathers seem to fol- 
low a biased erratic motion with the constraint of not 
going back in the direction they came from. Comparing 
the scaling of the diffusion associated with this peculiar 
behavior with the classical Brownian motion and biased 
random walk variants thereof would be an interesting av- 
enue for future exploration. 

In general, multi-peak breathers, or collections of 
breathers, also become mobile in the lattices with block 
capacitors. In all of our experimental results, the rele- 
vant motion is complex yet the relative distance between 
peaks remains fairly constant and does not fall below a 
minimum value (two edges away) (cf. Fig. [5]). This obser- 
vation suggests that two ILMs tend to repel each other 
upon close proximity. We also note, as in the square lat- 
tice, that although most transitions during motion hap- 
pen through the lattice edges, there are some transitions 
that happen along the long diagonal of the honeycomb 
cell (cf. Fig.[S]). We note that we did not observe transi- 
tions along the short diagonals of the honeycomb cell. 

A different scenario takes place when a smaller block 
capacitance, Cf = 15 nF, is chosen. In that case, 
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FIG. 4: (color online) Trajectory of a one peak moving 
breather on a square lattice for Vd = 2 V, / = 300 kHz and 
Cj — 1 fjF. Depicted are the positions of the center of mass 
of the breather in the (n, m) plane (left) and (n, m, t) space 
(right). The small (red) dots depict the experimental results 
and the large (blue) dots depict the corresponding numerical 
simulations. The trajectories have been unwrapped from the 
periodic lattice. Solid lines have been added for guidance. 




FIG. 5: (color online) Left: Experimental two peak moving 
breather profile on a honeycomb lattice with periodic bound- 
ary conditions for Vd = 2.5 V, / = 330 kHz, and C/ = 1 
fiF. Right: Trajectory of the corresponding centers of mass of 
the two interacting breathers (circles [red] and squares [blue]). 
Solid lines have been added for guidance. 



breathers become mobile but quickly collapse into a col- 
lective pattern resembling a planar (one-dimensional) 
wavefront, which rapidly transports energy coherently 
through the system, as shown in Fig. [6] It is intriguing 
that a spatially homogeneous driver can sustain such non- 
linear collective patterns characterized by energy trans- 
port. We should point out that the types of coherent 
patterns experimentally observable in the honeycomb lat- 
tice are apparently more complex than in the square lat- 
tice where only true planar wavefronts are seen. For the 
honeycomb lattice the equivalent planar waves, because 
of the geometry of the lattice, look more complex, and 
can be found at particular driving frequencies and am- 
plitudes (results not shown here) . Numerical simulations 
performed on larger square lattices (results not shown 
here) suggest that one-dimensional wavefront profiles are 
not robust — presumably due to transverse instabilities — 
and break down into patterns consisting of several ILMs, 
each one extending over a number of sites. A more de- 
tailed account of the phenomenology in larger lattices is 
currently under investigation and will be reported in a 
future publication. 

Conclusions & future challenges. In summary, we have 
generated two-dimensional discrete breathers in the set- 
ting of dampcd-driven electrical lattices. To our knowl- 




FIG. 6: Experimental (left) and numerical (right) moving 
front profile for a 6 x 6 square lattice. Vd = 2 V, / = 300 kHz 
and C f = 15 nF. 



edge, this is the first time that breathers have been exper- 
imentally stabilized in higher-dimensional discrete lat- 
tices by direct (and subharmonic) driving, with motion 
systematically induced to them and that these features 
have been examined not only in square but importantly 
also in honeycomb lattices. 

We have characterized the statics, stability and mobil- 
ity of these modes of self-localized energy, centered over 
a particular node of the lattice, and extending over a few 
lattice sites. The breathers have been found to persist 
indefinitely, and are either stationary or hopping in the 
lattice, depending on the precise make-up of the unit cell. 
We have observed that interacting breathers seem to re- 
pel each other when placed in close proximity which has 
the effect of (roughly) maintaining their pairwise distance 
in small periodic lattices. 

This work paves the way for numerous studies. These 
include the detailed characterization of breathers, their 
stability and mobility properties and associated potential 
barrier in long (and infinite) 2D chains; the characteriza- 
tion of impurities and their role in inducing mobility; the 
examination of spectral gap and higher gap (nonlinear) 
states, especially in honeycomb lattices; and the gener- 
alization of such chains even in fully 3D configurations. 
These issues of broader interest to other settings such as 
Josephson-junction ladders [l3| and granular crystals [l4[ 
will be considered in future publications. 
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